Inversion doublets of 3iV+iV cluster structure in excited states of 4 He 
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Excited states of 4 He are studied in four-body calculations with explicitly correlated Gaussian 
bases. All the levels below E yL =2& MeV are reproduced reasonably well using realistic potentials. An 
analysis is made to show how the state becomes a resonance but those having almost the same 
structure as this state in different spin-isospin channels are not observed as resonances. The role of 
tensor force is stressed with a particular attention to the level spacing between the two 0~ states. The 
calculation of spectroscopic amplitudes, nucleon decay widths, and spin-dipole transition strengths 
demonstrates that the state and the three lowest-lying negative-parity states with 0~ and 2~ 
have 3 H+p and 3 He+n cluster configurations, leading to the interpretation that these negative-parity 
states are the inversion doublet partners of the state. 

PACS numbers: 27.10.+h, 21.10. Jx, 21.45.-v, 21.60.De 



I. INTRODUCTION 

The competition of particle-hole and cluster excita- 
tions is one of the most interesting issues in the struc- 
ture of light nuclei. They emphasize different aspects of 
nuclear excitation modes, and often coexist in low-lying 
spectrum. Both the excitations are usually described in 
quite different languages, thus defying the reproduction 
of such a coexistence in a single scheme. In fact some 
intruder states have remained unresolved even in a large- 
space calculation based on realistic interactions. For ex- 
ample, the excitation energy of the so-called Hoyle state 
which is recognized to have large overlap with 3a con- 
figuration [l[ is predicted too high in the no-core shell 
model 0. According to the shell model, negative-parity 
states should appear first in the excited spectrum of 16 O, 
but they show up just above the first excited + state 
which is also understandable from 12 C+cv structure 

The 4 He nucleus is a lightest system offering a coex- 
isting spectrum. Its ground state is doubly magic and 
tightly bound, but its first excited state is not a nega- 
tive parity but + , similarly to 16 O. This state was first 
conjectured as a breathing mode, but an extensive study 
has confirmed it as a cluster state of 3N+N ( 3 H+p and 
3 He+n) configuration Accepting this interpretation 
for this state, we are led to the following questions. Since 
the 37V and N clusters having spin 1/2 and isospin 1/2 
move in a relative S wave, four states may appear which 
all have basically the same 3N+N configuration but dif- 
ferent J*T with 0+0, 1+0,0+1, 1+1. These states may 
be called quadruplets. The first question we set here is 
'Why do we actually observe only one of them, 0+0?' 

The second question is concerned with the concept 
of an inversion doublet which is known in molecular 
spectroscopy. For a system consisting of asymmetric 
molecules (clusters), one may expect a partner state 
of negative parity as in the ammonia molecule. These 
positive- and negative-parity pairs are called inversion 
doublets. On the analogy of the molecule we may 
ask a question 'What about the possibility of observing 



negative-parity partners in which the 3iV and N clusters 
move in a relative P wave?' The negative-parity partners 
would have J K =0~, 1~, and/or 2~, which result from the 
coupling of the spins of the two clusters and the relative 
orbital angular momentum between them. The centrifu- 
gal barrier for the P wave is more than 3 MeV at the 3N- 
N relative distance of 4fm, so that the expected partner 
states may appear in the region of the excitation energy 
_E x =21-23 MeV. In fact the 0~ and 2~ states are observed 
in this region. Traditionally, these states are considered 
s 3 p shell-model states, but could be better understood 
from the 3N+N configuration. The ATMS variational 
calculation seems to suggest this picture for the negative- 
parity states [5( but no discussion was made on the re- 
lationship between them and the 0+0 state. According 
to recent large-space shell-model calculations, these 0+0 
states show quite different convergence [f| : Because of its 
slow convergence, the 0^~0 state is attributed to a radial 
excitation. 

The purpose of this study is to answer the two ques- 
tions by performing four-body calculations with realistic 
potentials. Thus we are mainly interested in the three 
excited states, 0+0 (£ x =20.21 MeV), (TO (21.01 MeV), 
and 2~0 (21.84 MeV), but also consider other excited 
states which all have a width larger than 5 MeV. We 
will not invoke any model ansatz, that is, our calcula- 
tion is based on neither the shell model nor an RGM 
calculation 0, H[ which couples 3 H+p, 3 He+n, and d+d 
two-cluster channels, but treat four nucleons equally in 
an unconstrained configuration space. We will obtain 
the energies and wave functions of the excited states of 
He in a basis expansion method. The basis used here is 
square integrable, so that the excited states are obtained 
in a bound-state approximation. As we show later, this 
approximation works fairly well for predicting the three 
lowest-lying excited states, but it gives only a qualitative 
prediction for the other broad levels. 

SectionED gives a brief description of the basis functions 
used to solve the four-body problem. Section ITTTl presents 
the results of calculation together with some discussions. 
We show the energy spectrum of 4 He in Sec. IIII A[ dis- 
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cuss the problem about the quadruplets in Sec. IIII 51 and 
answer the question concerning the inversion doublets in 
Sec. IIII Cl Section llVl draws a conclusion of the present 
work. 



II. FORMULATION 



The Hamiltonian H for a system of two protons and 
two neutrons consists of the kinetic energy (T) and a 
nucleon-nucleon potential including the Coulomb poten- 
tial (Vcoui)- The center of mass kinetic energy is prop- 
erly subtracted. A three-body force is ignored as it has 
a small effect on the spectrum above the 3N+N thresh- 
old [|. We use the G3RS potential QjJ and the AV8' 
potential [ll[ as the two-nucleon interaction. Both of 
them contain central (V c ), tensor (Vt) and spin-orbit (Vt>) 
terms. The L 2 and L ■ S terms of the G3RS potential 
are ignored. The ground-state properties of d, 3 H, 3 He, 
and 4 He given by these potentials are similar to each 
other [12| . The tensor and spin-orbit forces of the AV8' 
potential are, however, stronger than those of the G3RS 
potential, while the central force of the AV8' potential is 
weaker than the one of the G3RS potential. 

A variational solution ^jMjTMt f° r the Schrodinger 
equation is obtained by taking a linear combination of 
many basis states, each of which has the following form 



T=0 and 1: 

,r {ls) 

0+ (00), (22); (11) 

1+ (01), (21), (22); (10), (11), (12), (32) 

0- (11); (22) 

1- (10), (11), (12), (32); (21), (22) 

2- (11), (12), (31), (32) ; (20), (21), (22), (42). 

Here the semicolon divides a natural parity set from an 
unnatural parity one. The values of L\, L2 in Eq. CD for 
a given L are chosen to be L, for the natural parity and 
L, 1 for the unnatural parity. Any basis functions with 
L 7r =0~ are not included in the present calculation. 

Each basis function differs in the choices of A, u\, and 
U2- The exponential part specified by A is called an ex- 
plicitly correlated Gaussian. An alternative expression 
for this part is given using the single-particle coordinate 
r, as 111 



exp 



2 ^ 

2<j 



V bu > 



(2) 



Specifying the elements of A using the 6 variables, 
(^12, ■ ■ • , ^34)7 is convenient for controlling the spa- 
tial extension of the system. 



III. RESULTS 



^(LS)JMjTMt 

= a{c-^ Ax [[y Ll (^)yL 2 (^)] L xs] JMj VTM T }, (1) 

where J^(r) = r e Yg(r) is a solid spherical harmonic. Here 
A is the antisymmetrizer, x a column vector whose ele- 
ments are three relative coordinates (cci, X2, X3), and A 
is a 3x3 positive-definite, symmetric matrix whose 6 in- 
dependent elements are variational parameters. The vec- 
tors iti and U2 each contain three elements determining 
the weightings of the relative coordinates and are used to 
specify the angular motion of the basis CD- The tilde " 
stands for the transpose of a column vector, and thus the 
inner product uix, which we call a global vector, is a vec- 
tor in three-dimensional coordinate space. However, the 
inner product x Ax denotes a scalar in three-dimensional 
space as it is defined by Xi-(Ax)i= J2i j AijXi-Xj. 

The global vector representation for the rotational mo- 
tion used in Eq. CD is found to be very useful. A reader is 
referred to Refs. [lllljl for 

more detail. The spin func- 
tion xsm s m Eq. CD is specified in a successive coupling, 
[[[||]si2 |]si23 |]sm s ! and all possible intermediate spins 
(S12, S123) are taken into account in the calculation. The 
isospin function t]tm t 1S Sl ^ so treated in exactly the same 
way as the spin function. We include the following (LS) 
values in Eq. CD to obtain the state with J* for both 



A. Energy spectrum 

The accuracy of our solution depends on the basis di- 
mension and the optimization of the variational parame- 
ters. The selection of the parameters is performed by the 
stochastic variational method [3, [H| ■ As all the states 
but the ground state are resonances, increasing the basis 
size unconditionally does not always lead to a solution 
we are seeking. Namely, if the variational parameters are 
allowed to reach very far in the spatial region, the en- 
ergy for the excited state would fall down to the 3 H+p 
threshold. 

Some details of the calculation are given below. The 
bij parameters are restricted to < fry < 8 fm for all the 
states but the 0~1 state. This choice covers the configu- 
ration space large enough to obtain accurate solutions for 
both the ground and first excited states 12]. Each ele- 
ment of Ui is allowed to take a value in the interval [— 1 , 1] 
under the constraint that its norm is unity, i.e. UiUi—l. 
Note that changing the normalization of Ui does not ac- 
tually alter the basis function CD except for its normal- 
ization. We found that using 600 basis states of the form 
given in Eq. CD (that is, 600 choices of parameters for 
A, ui, w 2 , L, S, S12, 5*123, T12, T123) enabled us to ob- 
tain converged solutions for both the ground and first 
excited states. See Fig. 1 of Ref. [l2j. Solutions for the 
other states are obtained in the basis dimension of 300. 
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FIG. 1: The energy convergence of the three lowest-lying 
negative-parity states of 4 He calculated using the G3RS po- 
tential. 



about 0.7 MeV with the experimental one. The theo- 
retical spectrum in the figure is drawn by shifting the 
calculated energies downward so as to adjust this dif- 
ference of the threshold energy. The theory reproduces 
the level sequence of the spectrum as a whole and espe- 
cially the excitation energies of the + 0, 0~0, and 2^0 
states very well. The levels above i? x =23MeV are pre- 
dicted to be slightly lower than experiment except for the 
0~1 level with AV8'. As their widths are all larger than 
5 MeV, this discrepancy may be allowable in the bound- 
state approximation for unbound states. Noteworthy is 
that the calculation predicts three states with + l, l + 0, 
and 1+1 around E x =23 MeV, as denoted by dashed lines. 
These states together with the + state are the quadru- 
plets relevant to the first question. As speculated, they 
show up in the present calculation, but no such states are 
observed experimentally. 



Figure [T] displays the energy convergence of the three 
lowest-lying negative-parity states with JT=0 _ 0, 2~0, 
and 2~1 as a function of basis dimension. The energies 
of these states are considerably stable for the increase of 
basis functions, though they do not have proper asymp- 
totic behavior characteristic of a resonance. The ener- 
gies of the other levels are not very accurately obtained. 
They have a large width, so that we think our calculation 
gives only an approximate energy. In particular we found 
that the energy of the 0~1 state, which has a width of 
about 8 MeV, was not as stable as the other states. We 
thus obtained its energy by restricting the range of fry as 
< bij < 6 fm. 

Figure [2] compares the spectrum of 4 He between theory 
and experiment. The calculated binding energy of 3 H is 
7.73 MeV for G3RS and 7.76 MeV for AV8' 12]. Thus 
the calculated 3 H+p threshold energy does not agree by 
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FIG. 2: Energy levels of the excited states of 4 He labeled with 
JT. Three of the quadruplets are drawn by dashed lines. 
The dotted lines indicate the 3 H+p and 3 He+n thresholds, 
respectively. Experimental values are taken from Ref. [16| . 



B. Quadruplets 

To resolve the first problem on the quadruplets, we 
begin by understanding why only the + state gets con- 
siderably lower than the other quadruplet members. As 
shown in Table[H all the members of the quadruplets con- 
sist of about 93% L=0 and 7% L—2 components. These 
values are almost equal to the corresponding components 
in 3 H and 3 He [12|, consistently with the conjecture that 
the quadruplets have 3N+N cluster structure with a rel- 
ative 5-wave motion. 

We list in Table |TT] the energy contents of the quadru- 
plet members. The matrix elements of the Hamiltonian 
H as well as its every piece are shown in a matrix form. 
The row and column labels of the 3x3 matrix correspond 
to the configuration space with L=0, 2, and 1. The off- 
diagonal elements in the lower triangular part of the ma- 
trix are added to the corresponding elements in the upper 
triangular part. We see that the key elements which give 
only the 0^0 state about 3 MeV larger binding energy 
than the other members are the kinetic energy as well as 
the tensor force. The kinetic energy contribution from 
the main channel listed in Table U is found to be about 
2 MeV smaller in the 0^0 state than in the other states. 
This is a consequence of the symmetry of the orbital part 
of the wave function as understood from Wigner's su- 
permultiplet theory [17]- The spin and isospin function 



TABLE I: Percentages of (LS) components of the quadru- 
plets calculated using the G3RS potential. Unnatural parity 
components are negligibly small. 





0+0 


0+1 


1+0 


1+1 


(00) 


93.0 


93.4 






(01) 






93.3 


93.4 


(21) 






3.0 


3.4 


(22) 


6.9 


6.6 


3.7 


3.1 
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TABLE II: Total energy of the quadruplets, given in MeV, and its decomposition into the contributions from the kinetic energy 
and the different potential pieces. The row and column of each 3x3 matrix correspond to the configuration space with L=0, 2, 
and 1. The off-diagonal elements in the lower triangular part of the matrix are added to the corresponding elements in the upper 
triangular part. For example, 4.58, —22.65 and —0.00 in the first row means the energy contribution of (L,L')=(0,0) channel, 
(0,2) and (2,0) channels and (0,1) and (1,0) channels, respectively. The contributions of L=3 are negligible and omitted. The 
G3RS potential is used. 
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4.58 


-22.65 


-0.00 


6.48 


-21.74 


-0.01 


6.30 


-21.67 


-0.00 


6.62 


-21.31 


-0.01 


(H) 




10.97 


-0.29 
0.14 




10.64 


-0.16 
0.08 




10.58 


-0.12 
0.06 




10.45 


-0.15 
0.08 




29.26 






31.09 






31.45 






31.73 






(T) 




10.99 


- 
0.15 




10.30 


- 
0.08 




10.31 


- 
0.06 




10.08 


- 
0.08 




-25.07 






-25.00 






-25.54 






-25.50 






(Vc) 




-1.56 


- 

—0.01 




-1.26 


- 

—0.00 




-1.28 


- 

—0.00 




-1.17 


- 

—0.00 




0.39 






0.39 






0.39 






0.40 






(V C oul) 




0.03 


0.00 




0.03 


0.00 




0.03 


0.00 




0.02 


0.00 






-22.65 






-21.75 
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-21.31 




(Vt) 




1.54 


-0.30 
0.01 




1.61 


-0.16 
0.00 




1.56 


-0.13 
0.00 




1.55 


-0.13 
0.00 








-0.00 






-0.01 






-0.00 






-0.01 


(Vb) 




-0.02 


0.00 
0.00 




-0.03 


0.00 
0.00 




-0.05 


0.00 
0.00 




-0.03 


0.00 
0.00 



of four nucleons contains more number of antisymmetric 
pairs in S=0, T—0 channel than in other ST channels, so 
that the orbital part of the 0^0 state is more symmetric 
with respect to the nucleon-exchange than the other or- 
bital functions. Furthermore, the 0^0 state gains about 
1 MeV energy compared to the others states, due to the 
tensor coupling between the L=0 and L=2 components. 

Now we discuss whether or not the quadruplets can 
be observed as resonances in 3 H+p and 3 He+n decay 
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FIG. 3: (Color online) SAs of the quadruplets for the S-wave 
3 He+n decay. The G3RS potential is used. 



channels. To this end we calculate a spectroscopic (or 
reduced width) amplitude (SA) defined as 



y{r) 



5(R 



Rr 



(N) Y t (R) 



JMj 



(3) 



Here R is the 3N-N relative distance vector, *i/2,i/2mi 
the normalized 3iV ground-state wave function, and 
^1/2,1/2— m t i s the nucleon spin-isospin function. They 
are coupled to the channel spin /. The label m-t dis- 
tinguishes either 3 H+p (mt=l/2) or 3 He+n (mt—— 1/2) 



channel. The 37V wave function ^1/21 



/2,l/2m, 



used here is 



obtained in the calculation using the basis of type (Q} 
with (LS)=(0±) and (2§) [12]. The component of L=l 
is very small (0.05%) and ignored unless otherwise. Fig- 
ure [3] displays the 3 He+n SA of the quadruplets. The 
orbital angular momentum between the clusters is set to 
£=0, and so / is equal to J. The 3 H+p SA is virtu- 
ally the same as the 3 He+n SA (except for the phase). 
The 0^0 state exhibits behavior quite different from the 
others: Its peak position is 3 fm, outside the 3A^ radius 
(~2.3fm). Moreover, the spectroscopic factor defined as 
J °° y 2 (r)r 2 dr is so large as 1.03. In a sharp contrast to 
the 0^0 state, the SAs of the other quadruplet mem- 
bers show nothing of resonant behavior: The peaks are 
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located extremely far outside the 3A radius, and the y 2 
value is small in the inner region. In passing we note that 
the SA of the ground state has a sharp contrast with that 
of the 0^ state: The peak appears near the origin and 
the amplitude is confined mostly in the 3 A radius. 

Since our variational solution is expected to be fairly 
accurate at least in the inner region, a decay width can 
be estimated with the formula of i?-matrix type: 

T N = 2P ( (kr)^-y 2 (r), (4) 



where k is the wave number given by k=y^2^iE /K 2 with 
the decay energy E, /i the reduced mass of the decaying 
particles, and Pi is the penetrability 

w= j ^)%m (5) 

n \2\ in \ i — in m for °He+n 

which is expressed in terms of either Coulomb wave func- 
tions or spherical Bessel functions. The decay width 
depends on the channel radius r, but its dependence is 
found to be mild: The T p value of the 0^~0 state is 0.69, 
0.74, 0.67 MeV at r=4, 5, 6fm, in good agreement with 
the empirical value of 0.50 MeV 16]. The above anal- 
yses all confirm that the 0^0 state has well-developed 
3N+N cluster structure, in accordance with the conclu- 
sion of Ref. [j]. Moreover, we conclude that none of the 
quadruplets except for the 0^0 state is a physically ob- 
servable resonance. This conclusion is consistent with 
the RGM phase-shift analysis which finds no resonance 
around 23 MeV excitation energy regionQ. 

As discussed above, all the quadruplet members but 
the 0^~0 state do not gain energy large enough to come 
down below the 3 He+n threshold. Only the 0^0 state 
shows up between the H+p and 3 He+n thresholds 
thanks to their Coulomb energy difference. The isospin 
conservation gives the 0jj~0 state an almost equal mixing 
of the open ( 3 H+p) and closed ( 3 He+n) channels. Both 
effects of the isospin conservation and the 3 H+p Coulomb 
barrier make the F p value of the 0^~0 state rather small. 
This state is thus a good example of a Feshbach reso- 
nance fisl ]. 

C. Negative- parity partners of the first excited + 

state 

Before coming to the inversion doublet issue, we first 
comment on the features of the negative-parity states. 
According to the shell model, the negative-parity states 
basically arise from the s^^ 2 p3/2 or s^ 2 p\ji particle-hole 

excitation, which predicts J^—O^, 1~, 1~, and 2~ for 
both T=0 and 1. However, a suitable combination of 
the two l - states with T—0 corresponds to the excita- 
tion of the center of mass, leaving only one 1 _ state with 
T=0. Seven negative-parity states observed experimen- 
tally below £? x =26MeV include three states with T=0 



TABLE III: Percentage of (LS) components of the negative- 
parity states calculated using the G3RS potential. The natu- 
ral and unnatural parity channels are separated by the line. 
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T=0 






T-- 


= 1 






0~ 


1" 


2- 


0~ 


17 

1 




2- 


(10) 




19.7 






51.0 


42.9 




(11) 


95.5 


74.2 


93.0 


96.9 


43.0 


53.1 


93.7 


(12) 




0.8 


0.3 




0.0 


0.6 


0.2 


(31) 






2.9 








2.8 


(32) 




3.4 


2.0 




4.3 


0.1 


1.7 


(20) 






0.0 








0.0 


(21) 




1.8 


0.5 




1.1 


1.5 


0.5 


(22) 


4.5 


0.2 


1.4 


3.1 


0.5 


1.8 


1.1 


(42) 






0.0 








0.0 



and four states with T=l, which is in agreement with 
the shell-model prediction. However, this agreement may 
not necessarily mean that the negative-parity states have 
shell-model like structure because the present four-body 
calculation also produces seven negative-parity states, as 
shown in Fig. [2] 

The level sequence is 0~, 2~, and l - in the order of 
increasing energy for T=0, while it is 2~, 1~, _ , and 1~ 
for T=l. Therefore the energy difference between the 
0"0 and 0~1 states becomes much larger than the one 
between the 1~0 and 1~1 states or between the 2~0 and 
2~1 states. It is interesting to clarify the mechanism of 
how this large energy difference is produced compared to 
the other negative-parity states with the same J w . Ta- 
ble [Hi] lists the percentage analysis of the seven negative- 
parity states according to their (LS) channels. The main 
component has L=l as expected from the shell model. 
The percentages are rather similar between the states 
with the same J™ but different T values. The similar- 
ity of the percentages is, however, not very clear in the 
1~ states because the values are fragmented into the two 
1~1 states. The main channel with L=l itself has a con- 
tribution from the tensor force but also gets a contribu- 
tion from the other channels through the tensor coupling. 
For example, the tensor force couples the natural parity 
channel (11) with the unnatural parity channel (22). 

Table llVl lists the energy content contributed from each 
piece of the Hamiltonian. Most striking is a different 
contribution of the tensor force. Compared to the 0~1 
state, the 0~0 state gains about 9 MeV energy from the 
tensor force, while the contribution of the central force 
to the energy gain is only about its half. The energy 
contents given by the AV8' potential are similar to those 
of the G3RS potential: The gain by the tensor force is 
even larger, about 12 MeV and the central force gives 
3 MeV difference. The tensor force is most attractive in 
the triplet even NN state, and it can be taken advantage 
of by having more number of antisymmetric AW pairs in 
the isospin space. The number of such antisymmetric 
pairs is counted from (r] TMT | J2i<, ( l " T i ' T j)/MVTM T ), 



6 



TABLE IV: Energy contents, given in MeV, of the negative-parity states. The G3RS potential is used. 

(TO 0~1 2~0 2~1 1~0 1~1 1J1 

(H) -6.40 -2.86 -5.78 -4.62 -3.69 -4.54 -1.84 

(T) 48.38 39.10 41.08 40.25 37.72 39.30 32.48 

(V c > -28.92 -24.79 -25.71 -25.82 -23.50 -25.14 -22.01 

(Vcoui) 0.48 0.44 0.42 0.43 0.40 0.43 0.42 

(K) -26.63 -17.75 -21.39 -19.30 -18.32 -19.13 -12.67 

(Vb) 0.29 0.14 -0.18 -0.18 0.006 0.005 -0.06 



which is [A(A + 2)-4T(T+l)]/8 for A-nucleon system. 
Thus the 0~0 state gains more attraction than the 0~1 
state through both the (ll)-(ll) diagonal and (ll)-(22) 
off-diagonal contributions [12j. If the unnatural parity 
basis were not included in the calculation, the (Vt) value 
of the 0~0 state would decrease to about half [12| and the 
0~0 state would lose significant energy. The role of the 
tensor force in lowering the energy of the 0~0 state was 
discussed many years ago . To be exact, the energy 

difference between the two states is actually a combined 
effect of the tensor, kinetic and central terms. 

Now we discuss the characteristics of the low-lying 
negative-parity states from the viewpoint of clustering. 
In Fig. 2] we display the 3 He+n SAs calculated from the 
three lowest-lying negative-parity states with 0~0, 2~0, 
and 2~1. The £ value for the 3 He-n relative motion is 
1, and the channel spin / is 1. As expected, each of 
the three curves shows behavior suggesting 3N+N clus- 
ter structure: The peaks are centered around 2 fm near 
the 3 AT surface, and the y 2 values are fairly large in the 
inner region. It is the centrifugal potential that makes 
the peak positions closer to the origin than that of the 
first excited + state. The 3N+N spectroscopic factors 
are considerably large: 0.58, 0.52, and 0.53 for the 0~0, 
2"0, and 2~1 states. We estimate the nucleon width. Its 
channel-radius dependence is again mild, so we choose 
r=5fm. The results for (T, r p /r), where T is the total 
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FIG. 4: (Color online) SAs of the three lowest-lying negative- 
parity states for the P-wave 3 He+n decay with 1=1. The 
G3RS potential is used. 



width in MeV, are (0.61, 0.72), (1.14, 0.58), and (1.85, 
0.53) for 0~0, 2~0, and 2~1, respectively. These val- 
ues are compared to those extracted from the P-matrix 
analysis [3, (0.84, 0.76), (2.01, 0.63), and (5.01, 0.53). 
The theory predicts the width of the 0~0 state very well, 
and gives about half of the width for the other states. 
Though the degree of clustering is somewhat reduced in 
the negative-parity states compared to the 0^~ case, the 
analysis of SA and decay width supports our conjecture 
that the 0~0 and 2~0 states as well as the 2~1 state con- 
stitute inversion-doublet partners of the first excited + 
state. The RGM phase-shift analysis of 3 H+p scatter- 
ings @ supports the P-wave resonance interpretation for 
these negative-parity states. The SAs of the 1~ states 
around P' x =24MeV show some degree of 3N+N clus- 
ter structure, though their amplitudes are considerably 
small compared to the 0~0 and 2~0 states in particular. 

An inversion doublet picture in nuclei was first pro- 
posed to understand the low-lying positive- and negative- 
parity rotation bands in 16 and 20 Ne from an a-core 
molecular structure [2(J. The appearance of positive- 
and negative-parity partners is a natural consequence of 
the underlying intrinsic structure dominated by the exis- 
tence of asymmetric clusters. We have shown that the 
three lowest-lying negative-parity states have a signif- 
icant component of 3N and N clusters whose relative 
motion is in P-wave. It is important to realize that this 
result is obtained in the calculation which assumes no 
cluster ansatz for the wave functions. A physical reason 
for the appearance of the inversion doublet partners is 
that they are located near the 3N+N threshold. 

Very unique in the inversion doublets in 4 He is that 
the 3N and N clusters have both J =1/2, and the chan- 
nel spin / is different in the doublets: It is for 0^0 and 1 
for 0~0, 2~0, and 2~1. The negative-parity partners with 
T=0 should thus be characterized by the transition of an 
isoscalar spin-dipole operator, 0\„=J2i = i [cjX^i - X4)] \ v , 
where X4 is the center of mass of 4 He. Note that r^ — X4 
is proportional to the distance vector between nucleon i 
and the center of mass of the other three nucleons. The 
transition strength to the 0^0 state, |(OjO||0 o ||O"O}| 2 , 
is 11.9 fm 2 , which is 6.9 times larger than that to the 
ground state. Moreover, the strength |(0^0||O ||0"0)| 2 
between the doublet partners occupies 58% of the "sum 
rule" J2 n I <0+0| |C> | |0~0) | 2 , where n takes all 600 eigen- 
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states with + 0. A similar enhancement occurs for the 
2-0 state as well. The value |(0+0||O 2 ||2-0)| 2 /5 is 
21.7 fm 2 , which is about 24 times larger than the one 
to the ground state, and it corresponds to 78% of the 
total sum E„l(O+O||0 2 ||2-O>| 2 /5. 

For the transition between the 2~1 and 0^0 states, 

an isovector spin-dipole operator, Ca^.io— Si=i [ <T i x ( T *i — 
x^jx^Ts^ must be considered. The transition strength 
(OjO|||02,i|||2 _ l)| 2 /T5 is 17.4fm 2 , which is 16 times 
larger than that to the ground state, where the triple 
bar 1 1 1 indicates that the reduced matrix element is taken 
in both the angular momentum and isospin spaces. This 
strength between the 2~1 and 0^0 states occupies 87% 
of the total strength |(0+0|||O 2 ,i|||2-l)| 2 /15. 

The high collectivity of the spin-dipole strength 
strongly indicates that the intrinsic structure of the 
negative-parity states, 0~0, 2~0, and 2~1, is similar to 
that of the first excited 0^0 state. 

IV. CONCLUSION 

A rich spectrum of 4 He comprising the coexisting levels 
has been reproduced in a single scheme without recourse 
to a specific model assumption. This has offered a good 
example of demonstrating the power of the global vector 
representation for the angular part of few-body systems. 



We have explained how only the 2 state is observed as a 
resonance among the quadruplets by examining the sym- 
metry property of the wave functions as well as the role 
of the tensor force. Analyzing the spectroscopic ampli- 
tudes, nuclcon decay widths and spin-dipole transition 
probabilities, we have confirmed that both the 0^ and 
negative-parity states with 0~0,2~0, and 2~1 are dom- 
inated by the 3N+N cluster structure and that these 
negative-parity states can be understood as the inversion 
doublet partners of the 0^" state in a unified way. Wc 
have shown that the tensor force plays a vital role to pro- 
duce the level spacing between the 0~ states with T=0 
and 1 through the coupling between the main channel 
with L=l and the unnatural-parity channel with L=2. 
A study of 16 in the scheme of 12 C+four nucleons will 
be interesting as its spectrum has some similarity to that 
of 4 He. 
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